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ABSTRACT 

We present an extinction map of a ~ 1 700 deg sq region that encloses the Ophiuchus, the Lupus, and the Pipe dark 
complexes using 42 million stars from the Two Micron All Sky Survey (2MASS) point sou rce catalog. The use of a 
robust and optimal near-infrared method to map dust column density (Nicer, described in Lombardi fc Alves]|2001 1 



allow us to detect extinction as low as Ak = 0.05 mag with a 1-a significance, and still to have a resolution of 3 arcmin 
on our map. We also present a novel, statistically sound method to characterize the small-scale inhomogeneities in 
molecular clouds. Finally, we investigate the cloud structure function, and show that significant deviations from the 
results predicted by turbulent models are observed. 

Key words. ISM: clouds, dust, extinction, ISM: structure, ISM: individual objects: Pipe molecular complex, Methods: 
data analysis 



1. Introduction 

In this paper we present an extinction map of the 
Ophiuchus and Lupus complexes covering ~ 1 700 deg sq, 
computed by applying an optimized multi-band tech- 
nique dubbed Near-Infrar ed Color Excess Revisited (Nicer 
Lombardi & Alves 200 1[ hereafter Paper I) to 42 million 
JHK photometric measurements of stars from the Two 



dition, [Knude fc Nielsen (20011 suggested that Lupus 2 is 
physically disconnected from the other Lupus subclouds, 



and is located at a distance of 360 pc; similarly, Bertout 



Micron All Sky Survey (2MASS; |Kleinmann et al.| [1994). d 



et al.| ( 1999) suggested Lupus 3 to be disconnected from the 
rest of the Lupus complex. However, surprisingly no veloc- 
ity difference is found among the various Lupus subclouds 
in radio observations. Recently, we used a novel maximum- 
likelihood techniqu e based on Hipparcos data obtaining 
(155 ± 8) pc (|Lombardi et al. 



2008); this value will 



This paper is the second of a series where we apply the be adopted in ^ pap er as the actual dis tance of Lupus 



Nicer method to point sources from the 2MASS database 
The main aim of this coordinated study is to investigate in 
detail the large-scale structure of molecular clouds. In ad- 
dition, the use of a uniform dataset and of a consistent and 
well tested pipeline allows us to characterize many proper- 
ties of molecular clouds, such as their reddening law, and 
to identify cloud-to-cloud variations in such properties. The 
region considered here encloses in addition the Pipe nebula, 
consi dered in the first paper of this series ( Lombardi et al. 
[20061 hereafter Paper II), and allows us to study in detail 



the relationships among these molecular clouds. 

Lupus is a well studied complex, composed of several 
subclouds showing different modes of star formation: Lupus 
1 shows isolated star formation, while cluster formation is 



observed in Lupus 3 


Schwartz 


1977 


Hughes et al. 


1994 


Nakajima et al. 


2000 


; no evident star formation activity 



plex is still highl y con troversial. The conventional value 



150 pc (|Krauuter 19911, has been challe nged by Knude & 
Hog ( |1998p , who reported 100 pc, and by Wichmann et a 



(1998J), 190 pc, both using Hipparcos observations. In ad 
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Due to its large area on the sky, millimeter observations 
have mostly covered limited regions of this cloud. The first 
ext ensive cloud su r vey in 12 CO (J = 1 — > 0) was made 



by Murphy et al. (19861 over ~ 170 sq deg with an ef- 
fective reso lution^ of 30 arcmin. More recently, [Tachihara 
et al. ( 2001 1 have performed a complete survey of the Lupus 
complex^also in 12 CO, J = 1 — * 0) using the NANTEN 
sub-millimeter telescope. Their data covered 550 sq deg on 
a grid of 8 arc min, and smaller area s on a 4 arcmin grid. 
Finally, recently |Teixeira et al. (2005) studied the structure 
of Lupus 3 in the NIR. 

The p Ophiuchi molecular cloud is a filamentary com- 
plex located at the edge of the Upper Scorpius subgroup 
in the Sco-Cen OB association. It is composed of a series 
of dark clouds that extend eastward from cores of dense 

With an estimated distance 



molecular gas (deGeus 



19921 



of (119±6) pc (Lombardi et al. 2008), it is one of the nearest 
star- forming regions, and has thus been the target of nu- 
merous investigations. The main cloud, L1688, is situated 
about one degree south of the star p Ophiuchi and has been 
studied in the near and far infrared, in the millimeter con- 
tinuum, as well as in the X-ray and radio continuum, and 
hosts an embedded infrared cluster with ~ 200 young stel- 
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lar objects (YSO; e.g., 

Comeron et al 



et al. 


2001 


L168b 


is kno 



Greene & Youn 



ig iyyzjjaontemps 
1'achinara et al ||2000p 
z central molecular cor< 



1992 Bontemps 



L1688 is known to host a 1 pc x 2 pc central molecular core 
with visual estimated v isual extinction exceeding 50 mag 



(Wilking & Lada 
an area of rough 



19831. This "main cloud" of Oph covers 
y 480 sq arcmin and has been dissected 
into about a dozen cloud components. The entire star for- 
mation complex extends over s everal degrees on th e sky, 



Motte et alT] ( |1998[ ) have 



containing a few major clouds. 

conducted an extensive 1.3 mm continuum mapping of the 
centra l region using the IRAM 30-m telescope. |Stanke et al.| 
(20061 presents SIMBA 1.2 mm observations of the cloud 
covering an area of 460 arcmin 2 with a resolution of 24"; 



Johnstone et al.| ( |2000| | have published a SCUBA 850 jum 
thermal emission map of the p Ophiuchi cloud. 

This paper is organized as follows. In Sect. [2] we briefly 
describe the technique used to map the dust and we present 
the main results obtained. A statistical analysis of our re- 
sults and a discussion of the bias introduced by foreground 
stars and unresolved substructures is presented in Sect. [3] 
Section |4] is devoted to the mass estimate of the cloud com- 
plex. Finally, we summarize the results obtained in this 
paper in Sect. [5] 



2. Nicer extinction map 

We carried out the data analysis using the Nicer method 
described in Paper I. Near infrared J (1.25 /im), H 
(1.65 /jm), and K s band (2.17 /im) magnitudes of stars in 
a large region of the sky which includes the Ophiuchus and 
the Lupus dark clouds were taken from the Tw o Micron All 
Sky Surve^f] (2MASS; iKleinmann et ai]|l994[ ). In particu- 
lar, we selected all 2MASS reliable point sources within the 
boundaries 



-32° < I < +16° 



-2° < b < +40° 



(1) 



This area is 1 672 square degrees and contains approxi- 
mately 42 million point sources from the 2MASS catalog. 
Note that since the region spans several degrees in galac- 
tic latitude, the local density of background stars changes 
significantly in the field; as a result, since we used a fixed 
size for the smoothing of the extinction map (see below), 
the noise of the map increases as b increases. 

When selecting stars from the 2MASS point source cata- 
log we discarded possible spurious detections (e.g., objects 
likely to be associated to solar system minor bodies, ex- 
tended sources, or artifacts). We then generated a prelim- 
inary extinction map which, as described in Paper I, was 
mainly used as a first check of the parameters adopted, to 
select a control region on the field, and to obtain the photo- 
metric parameters to be used in the final map (see Fig. pi . 
We identified a large region that is apparently affected by 
only a negligible extinction (see below) , and used the colors 
of stars in this control field as reference ones. 

Using the information provided by the control field, 
we generated a second map, which is thus "calibrated" 
(i.e., provides already, for each position in our field, a reli- 
able estimate of the colu mn density). Note that the recent 
Indebetouw et al. ( 2005 1 2MASS reddening law was used 
here. We then considered the color-color diagram for the 




-1 



Fig. 1. Color-color diagram of the stars in the whole field, 
as a density plot. The contours are logarithmically spaced, 
i.e. each contour represents a density ten times larger than 
the enclosing contour; the outer contour detects single stars 
and clearly shows a bifurcation at large color-excesses. The 
dashed lines identify the regions in the color space defined in 
Eqs. ^ and d3l, as indicated by the corresponding letters. 
Only stars with accurate photometry in all bands (max- 
imum l-cr errors allowed 0.1 mag) have been included in 
this plot. 
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Fig. 2. Spatial distribution of the samples of sources as de- 
fined by Eqs. p]) and Sample A is shown as filled circles, 
while sample £> is shown as crosses (see also Fig. [TJ . Sample 
A appears to be strongly clustered in high-column density 
regions of the cloud, and is thus interpreted as genuine red- 
dened stars; sample B seems not to be associated with the 
cloud, and is instead preferentially located at low galactic 
latitudes and in the Galactic bulge. 




See http : //www . ipac . caltech. edu/2mass/. 
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Fig. 6. The Nicer extinction map of the Ophiuchus, Lupus, and Pipe complexes. The resolution is FWHM = 3 arcmin. 
The various dashed boxes mark the regions shown in greater detail in Figs. [9] and 10 The overplotted contour is at 
Ak = 0.3 mag. The Pipe nebula occupies the region around (I, b) — (0°, 5°). 
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Fig. 3. The histogram of the K band magnitude for the 
two star subsets A and B of Eqs. ([2| and 



stars in the catalog to check for possible signs of anoma- 



lous star colors. The result, presented in Fig. [T] shows a 
bifurcation for H — K > 1 mag. 

As discussed in detail in Paper II, the bifurcation is 
likely to be created by Asymptotic Giant Branch (AGB) 
stars. To further investigate this point, we considered the 
two samples in the color-color diagram defined as 



A = {IA(H -K) + 0.5 mag < (J - H) 
and H — K > 1 mag} , 

B = {IA(H -K) + 0.5 mag > (J - H) 
and H — K > 1 mag} . 



(2) 
(3) 



An analysis of the spatial distribution of these two sam- 
ples (Fig. [2l reveals that, as expected, sample A is associ- 
ated with the densest regions of the molecular cloud, while 
sample B is distributed on the whole field with a strong 
preference for low galactic latitude regions. 

The nature of the two stellar populations in samples A 
and B is further clarified by the histogram of their K band 
magnitudes, shown in Fig. [3] As expected, sample A shows 
a broad distribution, which can be essentially described as 
a simple power-law luminosity function up to K ~ 12 mag; 
note that the completeness limit of our sample is signifi- 
cantly smaller than the typical 2MASS completeness in the 
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Fig. 7. The star density for the extinction map of Fig. ^ note how evident is the galactic bulge here. The faint white 
structures correspond to the position of dense cores, where the star density decreases because of the cloud extinction. 
The white "hole" at the center-bottom of the field is Antares. The black spots are open and globular clusters. 



K band (14.3 mag at 99% completeness) because of the 
stricter selection operated here (small photometric errors 
in all bands) and because most sample A stars come from 
low galactic latitude regions (where the increased density of 
stars significantly reduces the completeness of the 2MASS). 
In contrast to sample A stars, sample B stars show a well 
defined distribution, with a pronounced and relatively nar- 
row peak at K ~ 7 mag. This strongly suggests that we are 
looking at a homogeneous population of sources located at 
essentially the same distance. 

The lack of correlation between the dust reddening and 
the stars of sample B can be also investigated by consid- 
ering the extinction-corrected color-color diagram shown 
m Fig. |4] This plot was obtained by estimating, for each 
star, its "intrinsic" colors, i.e. the extinction corrected col- 
ors from the extinction at the star's location as provided 
by the Nicer map. In other words, we computed 

J intr = J — (Aj/A K )A K , (4) 

ff intr = H - (A H /A K )A K , (5) 

K intr =K-A K , (6) 

where Ak is the Nicer estimated extinction in the direc- 
tion of the star from the angularly close objects. By com- 



paring Fig. [4] with Fig. [T] we see that the upper branch, 
i.e. sample A, is strongly depressed, while the lower branch 
(sample B) is largely unaffected. Note that the residual 
stars appearing in the upper branch are likely to be the ef- 
fect of an inaccurate extinction correction due to small-scale 
inhomogcncities not captured by our analysis; similarly, the 
tail at negative colors is due to "over-corrected" stars (for 
example foreground stars observed in projection to a cloud). 
Because of the much shrunken distribution of the upper- 
branch, Fig. [4] is particularly useful to better identify stars 
belonging to the lower branch, and was thus used to define 
two further subsets: 

Bi = {\(J- H) iatr - IA(H - K) intI + 0.5| < 0.9 

and (J - fQintr > -4.6(# - K) iatl + 4} . (7) 

B 2 = {\(J- H) iatr - 1A(H - K) intl + 0.5| < 0.9 

and (J - H) intI > -2A(H - K) intI + 2.5} . (8) 

These two subsets, marked in Fig. [4j correspond to the ar- 
eas in the color-color plot where a contamination by sam- 
ple B stars is possible (i?i), or highly likely (B 2 ). We find 
approximately 996 000 B\ stars and 360 000 B 2 ones, cor- 
responding to 2.3% and to 0.88% respectively. 
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Fig. 8. The statistical error a{Ax) on the measured column density of Fig. Kil Note that correlation on the errors on a 
scale of FWHM = 3 arcmin is expected. The increase of the error with the galactic latitude is evident. 



In summary, all the evidence found supports the iden- 
tification of the "lower branch" with evolved intermedi- 
ate mass stars (about 1 to 7 M Q ), and likely with the 
Asymptotic Giant Branch (AGB) at about the distance to 
the Galactic center (see also Paper II). Since the "lower 
branch" stars seem to be unrelated to the molecular cloud, 
and since their colors would be interpreted by the Nicer 
algorithm as a sign of extinction, in principle our results 
could be biased toward a higher extinction, especially at 
low galactic latitude regions. In practice, we argue that the 
bias introduced by "lower branch" stars is negligible given 
the low density of these stars. 



Nevertheless, and in order to avoid any source of bias, 
although small, we excluded from the 2MASS catalogs all 
stars located in the B 2 color-space region, and performed 
the whole analysis described in this paper using this re- 
duced subset of stars. We stress that if we had performed a 
cut in the observed colors, we would have introduced a new 
bias in the deduced column density; instead, a selection in 
the intrinsic colors does not bias the final results. As an 
example, Fig. [5] shows the color-color diagram for the new 
set of stars: note that the "lower branch" disappears com- 
pletely in this plot, a further confirmation that our selection 
is effective in removing this population of stars. 



We then run again the whole Nicer pipeline on the re- 
fined catalog. After (re)evaluating the statistical properties 
of stars in the control field, we constructed the final map, 
shown in Fig. [6] We recall that in Nicer the final map 
can be generate d using different smoothing techniques (see 
Lombardi 2002 for a discussion on the characteristics and 
merits of various interpolators) . As pointed out in Paper I, 
generally these techniques produce comparable results, and 
thus we focused here to the simple moving weight average: 



A K (0) 



En=! 



(n) 
K 



(9) 



where Ak{0) is the extinction at the angular position 9, 

A$ is the extinction of the n-th star, and W ijL \6) is the 
weight for the n-th star for the pixel at the location 9: 



w {n \e) = 



w(e-e^) 

Var(i^) 



(10) 



Hence, the weight for the n-th star is composed by two fac- 
tors: (i) W(9 — 0( n '), i.e. a function of the angular distance 
between the star and the point 6 where the extinction has 

ljP), the inverse of the 



to be interpolated, and (ii) Var(A^ 
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Fig. 9. A zoom of Fig. [^showing the central region of the Ophiuchus cloud. The box encloses the core of cloud, for which 
a few separate studies are presented in this paper. 



inferred variance on the estimate of Ak from the star. The 
first factor, parametrized by the weight function W, is here 
taken to be a Gaussian. 

The map of Fig. [6] was generated on a grid of approxi- 
mately 2 000 x 1 600 points, with scale 90 arcsec per pixel, 
and with Gaussian smoothing characterized by FWHM = 
3 arcmin. Note that in the weighted average of Eq. Q we 
also introduced an iterative a-clipping at 3-cr error (see 
Paper I). The average, effective density of stars is ~ 6.5 
stars per pixel, but as noted above this value changes sig- 
nificantly on the field with the galactic latitude (see Fig. 7]); 
this density guarantees an average (1-c) error on Ak be- 
low 0.03 magnitudes. The largest extinction was measured 
close to p Ophiuchi, where Ak — 2.89 magnitudes. The ex- 
pected error on the extinction Ak is shown in Fig. [8] and 
was evaluated from the relation (see Paper II) 



4, 



(0) 



En=l[W {n H0)] ^v(A 



K ) 



(11) 



As expected, we observe a significant gradient along the 
galactic latitude. Other variations in the expected errors 
can be associated to bright stars (they produce the charac- 
teristic cross-shaped patterns), to globular clusters (white 
dots in Fig.|8|, and to the cloud itself (dark areas). Because 



of the relatively large variations on the noise of the extinc- 
tion map, clearly a detailed analysis of Fig. [6]should be car- 
ried out using in addition the noise map of Fig. [8} Figures [9] 
and 10 shows in greater detail the absorption maps we ob- 
tain for the Ophiuchus and Lupus complexes, and allow us 
appreciate better the details that we can obtain by applying 
the Nicer method to the quality of 2MASS data. 



3. Statistical analysis 

3.1. Reddening law 

The large number of stars in the field allowed us to accu- 
rately check the reddening law used throughout this paper. 
To the purpose, we partitioned all stars into different bins 

corresponding to the individual original Ajp measurements 
(we used a bin size of 0.02 mag). Then, we evaluated the 
average NIR colors in each group of stars in the same bin 
and the corresponding statistical uncertainties (estimated 
from on the photometry errors of the 2MASS catalog). The 
results obtained are shown in Fig. 11 together with the 
normal infrared reddening law in the 2MASS photometric 



system ( jlndebetouw et al.|2005 1 . This plot shows that there 
are no significant deviations from the normal reddening law 
over the whole range of extinctions investigated here. 
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Fig. 10. A zoom of Fig. [^showing the of Lupus cloud complex (in particular, all clouds from Lupus 1 to Lupus 9). 



3.2. Foreground star contamination 

If a fraction / of observed stars is foreground to a dark 
cloud, the measured column density is underestimated by 



a factor (1 — /) . This effect is normally negligible on the out- 
skirts of nearby molecular clouds, where the fraction / of 
foreground stars is typically of the order ~ 1%, but unfor- 
tunately can have significant effects on very dense regions, 
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-1 o 1 2 3 4 0.475 represent, respectively, the average of H — K and of 

(H-K) iat J — H colors in magnitudes for the control field). Error 

bars are uncertainties evaluated from the photometric er- 
Fig. 4. The extinction-corrected color-color diagrarm The rors of the 2 MASS catal og. The solid line shows the normal 
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Fig. 5. The color-color diagram for the selected stars in the 
field, after the removal of the set B2 of Eq. (j8| . A compari- 
son with Fig. [T] shows that we were able to virtually remove 
all significant contamination from spurious reddening. 



where because of a selection effect / increases significantly. 
In addition, many dense cores host young stellar objects: 
these stars, if moderately embedded, show only a fraction of 
the true, total column density of the cloud. As a result, the 
extinction in the direction of dense regions can be severely 
underestimated. For example, the core of p Ophiuchi shows 
an apparent "hole" in absorption (see Fig. 12), but a corn- 



Fig. 12. A zoom of the extinction map near p Ophiuchi, 
with overprinted the star density contours. Note the white 
"hole" close to the center of this core, which is due to the 
combined effect of very high extinction values and the pres- 
ence of embedded stars in this star forming core. 



parison with the local density of stars shows that the hole 
is the result of young stars moderately embedded on this 
active star forming region (this core host ~ 200 young stel- 
lar objects) and of foreground star contamination (see also 
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Fig. 13. The probability distribution of star pixel extinc- 
tions for the Ophiuchus map shown in Fig. [9] the gray, 
smooth curve represents the best-fit with a log-normal dis- 
tribution. 




Fig. 14. The probability distribution of star pixel extinc- 
tions for the subset b > 6.5°) of the Lupus region of Fig. 10 
the gray, smooth curve represents the best-fit with a log- 
normal distribution. 



Barsony et al. ( 1997 1 for a discussion on the effects of fore- 
ground stars in the p Ophiuchi core). 

In order to evaluate quantitatively the fraction / of fore- 
ground stars, we selected high-extinction regions character- 
ized by Ak > 1 mag. We then checked all stars in these 
regions that show "anomalous" extinction, i.e. stars whose 
column densities differ by more than 3-<r with respect to the 
field. A total of 1 268 stars met this criteria, but only 784 
of them show measured column densities compatible with 
no extinction. Hence, since the total area of regions with 
Ak > 1 mag is about 1.4 deg 2 , we estimate that on aver- 
age on the field only a fraction ~ 2% of stars is foreground. 
As a result, we can safely ignore the effect of foreground 
stars except on the higher extinction regions which, how- 
ever, are not the focus of this paper. 

3.3. Column density probability distribution 
Theoretical studies (e.g. Vazquez-Semadeni||1994| |Padoan 



et al.|1997b| Passot fc Vazquez- Scmadcni 1998; S calo et a 
1998 1 have shown that the probability distribution for 
the volume density in molecular clouds is log-normal for 
isothermal flows (i.e., when the polytropic index 7 = 
dlnP/dlnp = 1), while it tends to develop a power-law 



tail at high (respectively, low) densities for 7 < 1 (7 > 1). 
In reality, observations can only probe the probability dis- 
tribution for the column density, i.e. the volume density 
integrated along the line o f sight. As discussed by Vazquez- 1 
Semadeni & Garcia (2001 1, this quantity can behave differ- 
ently depending on the ratio r\ of the column density to a 
suitable defined "decorrelation length." In particular, if r\ 
is large, so that there are many independent "events" along 
the line of sight, the probability distribution for the column 
density is expected to be normal, a re sult in agreement w ith 



the central limit theorem (see, e.g., Eadie et al. 19711; if, 



instead, 77 is small, the probability distribution for the col- 
umn density is expected to be similar to the one for the 
volume density, i.e. log-normal in the isothermal case. In 
practice, numerical simulations show that the convergence 
to the normal distribution is extremely slow and that tran- 
sition between the two cases, in typi cal cases, is for 77 ~ 10 4 
( Vazquez-Semadeni fc Garcia||2001 1. 

Figures |13| and |14| report the probability distributions 
of column densities observed in the Ophiuchus and Lupus 
fields, i.e. the relative probability of column density mea- 
surements for each pixel of Figs. [9] and [To] We tried to fit the 
column density histograms with log-normal distributions of 
the form 



h(A K ) = 



A K - A Q 



exp 



(ln(A K -A Q )-lnA 1 
2(lna) 2 



2 1 



(12) 



For the Ophiuchus complex we obtained a satisfactory fit 
on the whole region; for Lupus, the fit is extremely good, 
especially when excluding the low-galactic latitude regions 
(b < 6.5°), which are likely to be contaminated by different 
cloud complexes. Note, by comparison, that as described 
in Paper I, for the Pipe nebula a good fit to the column 
density distribution requires four normal distribution. This 
result might indicate that the Pipe nebula is the result of 
the superposition of different components, each of which is 
likely to be extended on the line of sight (so that the central 
limit theorem can be applied to the total measured column 
density) . 

3.4. Small-scale inhomogeneities 

It has been long recognized ( Lada et al.|1994 ) that the local 
dispersion of extinction measurements increases with the 
column density: in other words, for a fixed (small) patch of 
the sky, the scatter of the individual star estimates of Ak 
increases as the average of Ak increases. The scatter in 
the Ak estimates is typically associated with the intrinsic 
scatter in the NIR star colors and to the effect of small- 
scale inhomogeneities in the cloud projected density. These 
inhomogeneities are important for a number of factors: 



Center Ao Scale A\ Dispersion Peak a 
Ophiuchus -0.081 Ol51 0681 946 

Lupus (b > 6.5°) -0.135 0.153 1.506 651 

Table 1. The best-fit parameters of the four Gaussian func- 
tions used to fi t th e col um n density pro bability distribution 
shown in Figs. 13 and 14 (see Eq. (12 1 for the meaning of 



the various quantities) 
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In the simplest interpretation, they indicates differen- 
tial extinction or strong gradients in the column den- 
sity (with typical scales smaller than the resolution of 
the extinction map). 

If present at different scales, they can be i nterpreted 
as the effects of turbulent motio ns (see, e.g. 
[Baily|[l994l |Padoan et al~|l997a| ). 
In presence of significant inhomogeneities at scales 
smaller than the resolution of the extinction map, 
Nicer (as well as other color-excess methods) is ex- 
pected to be biased toward low extinction. This hap- 
pens because the background stars will no longer be 
randomly distributed in the patch of the sky used to 
estimate the local extinction value, but will b e preferen- 
tially detected in low-extinction regions (see Lombardi 



2005 for a more detailed discussion on this point) 



In order to better quantify the effect of inhomogeneities 
on small scales, consider the quantity (cf. Paper II) 



(13) 



Note that <r^ is defined in a different way with respect 



to ct^ of Eq. (11). Let us now fix a given direction in the 



sky 9, and let us consider the process of measuring the col- 
umn density Ak(0) there. This quantity is evaluated using 
Eq. ([9]), i.e. Ak{0) is a weighted mean of the estimated 

column densities {A^ '} of the stars observed close to the 
direction 0. The estimated star column density for the n-th 
star, in turn, can be written as 



= A K 



A (n) 



An) 



(14) 



In this equation we split the three contributions to the 
measured column density: Ak, the "average" extinction 
in the patch of the sky considered (see below); AW = 
Ak(@^) ~ Ak, the local difference from the average ex- 
tinction, which includes both random inhomogeneities (e.g., 
due to turbulence) and unresolved structures (e.g., due to 
steep gradients in the extinction); and p n ', the photometric 
error on the measured extinction of the n-th star. 
Because of the presence of the photometric error 



in Eq. (14 1, we cannot use directly the estimator a - of 

-1 Ak 



Eq. (13 1 as a measure of small-scale inhomogeneities, and 



a more detailed analysis is needed. In the following we will 
consider ensemble averages with respect to the photometric 
errors e'"'; in other words, we will evaluate the mean values 
and variances of some relevant quantities by taking e^- n ' as 
independent random variables with the properties 



(eW) = , ((e< n) ) 2 ) = Var(i^) = V™ 



(15) 



where the last e qua lity is a mere definition. In the decom- 
position of Eq. ( 14 1 we can freely choose the value of the 



"average" extinction Ax (a change of Ak can be counter- 
balanced by a change of all A^™'); we now fix this quantity 
to 



A, 



ZLiW^A K {eW) 



(16) 



This is a convenient and natural choice since in this case 
we have 



0, 



(17) 



and (Ak(0)) — Ak- Note also that in this case we have 



A K {6)-A K 



E„=i W^e^ 



e{0) 



(18) 



The quantity e(8) can be interpreted as the "error" on the 
extinction map: 



(e(0)) = 0, 



<[^)] 2 > = 



(19) 
(20) 



where in the last step we used the definition (111 



We can now fina lly consider the ensemble average of the 
quantity of Eq. ( 13 I: 



SLi^ (n) ([AW + e W- e (0)] a 

\ a A K ) - yN w(n) 



Using the relation 

c (n) L-im " 



(e<">e(0)) 



(>") 



we can expand Eq. (21 1 into 
£„V^")(A(")) 2 



<*L> = 



v I ~ iu( m ) 



(21) 



(22) 



- 2 



E„^ (n) E„^ (n) 
E»(w (n) ) 2 v< B > 

(E„^(«)) 2 

E„ M/( " ) ( A( " ) ) 2 J2n W( - n)V( - n) 



4 K W 



(23) 



The first term in the r.h.s. of this equation represents a mea- 
sure of small scale inhomogeneities. Since all other quanti- 
ties appearing in Eq. ( 23 1 can be evaluated from the data, 



it makes sense to define the map 

J2 n w {n) v {n) 



E„^(") 



(») 



E„^ ( " } ( A 
E„^(«) 



(24) 



and to interpret it as a "variance" of small scale inhomo- 
geneities. Note finally that if the weights |W( n H are cho- 
sen according to Eq. ( 10 1, then in the numerator of the last 



term of this equation we can use W^V^ = W(0 — 

We evaluated the A 2 map for the whole field in order 
to identify regions with large small-scale inh omo gen eities. 
The outcome of this analysis (see Figs. 15 and 16 1 sup- 
port some results recently obtained for the Pipe nebula 
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Fig. 15. The A 2 map of E q. p4| on the same region around 
p Ophiuchi shown in Fig. |12[ with overplotted star density 
contours. Note the significant increase observed in A 2 close 
to the central parts of this cloud. 
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Fig. 16. The A map of Eq. (24) on a region close to 



Lupus 1, in logarithmic scale, with over plot ted start density 
contours (using the same levels as Fig. 15 1. 




A K (mag) 

Fig. 17. The distribution of the A 2 map as a function of 
the local extinction Ak for the Ophiuchus, in logarithmic 
grey scale. The dashed line shows the average values of A 2 
in bins of 0.05 mag in Ak- Note the rapid increase of A 2 
for Ak > 0.8 mag. As a comparison, the average variance 
Var(A^) on the estimate of Ak from a single star is ap- 
proximately 0.033 mag 2 . 




Fig. 18. Same as Fig. 17 for Lupus. Note that the effect of 
substructures is evident at significantly lower column den- 
sities. 



(Lombardi 2008[ ): inhomogeneities are mostly present in 
high column density regions, while at low extinctions (ap- 
proximately below Ak < 0.5 mag) substructures are ei- 
ther on scales large enough to be detected at our resolution 
(3 arcmin), or are negligible. The range of values spanned 
by the A 2 map is quite impressive: it reaches (1.5 mag) 2 in 
the p Ophiuchi core, indicating an extraordinary amount of 
small-scale substructures present there. The Lupus region 
shows less pronounced small-scale inhomogeneities, with a 
maximum A 2 ~ (0.66 mag) found in the Lupus 3 core (not 
shown here). The Lupus A 2 map at low galactic latitudes 
shows significant amount of substructures, possibly due to 
the presence of other molecular clouds (see below). 

Figures 17 and 18 shows the average A 2 as a function of 
the local extinction Ak for the Ophiuchus and the Lupus 
clouds. As shown by these plot, substructures start to play 
a significant role at relatively large column densities and 
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are negligible at low Ak- In particular, the dashed lines in 
these plots, representing the average value of A 2 in bins of 
0.05 mag in Ak, should be compared with the average vari- 
ance VarfAj^ ) on the estimate of Ak from a single star, 



17 



and 



18 



which is approximately 0.033 mag 2 . From Figs 
we thus see that local inhomogeneities start to be the preva- 
lent source of errors in extinction maps for Ak > 0.8 mag 
for the Ophiuchus cloud. For the Lupus complex, appar- 
ently this happens before, but a more detailed analysis 
shows that this result is mostly due to the rapid increase 
of A 2 at low galactic latitudes (which, in turn, is probably 
due to other intervening clouds). 

Diagrams such as the ones presented in Figs. [T7| and [T8| 
are invaluable to understand the small-scale properties of 
molecular clouds, but clearly they can only be interpreted 
with a detailed model. A thorough analysis of these results 
goes beyond the scope of this paper, and in any case is 
hampered by the relatively coarse resolution achievable us- 
ing 2MASS data. Still these data allow us to confirm that 
the observed scatter in column density is due to small-scale 
substructure and not to other effects such as foreground 
star contamination. The latter, typically, show up in the 
Ax-A 2 diagram as separate trails with parabolic shape that 
divert from the main A 2 = locus of points (cf. Fig. 9 in 
where however a » r ~ V A 2 is plotted, 



Lada et al. 



1994 



so that the parabolic shape becomes a line). This pattern 
typically is well discernible, and in absence of other sources 
of inhomogeneities produces a clear signature that is not 
compatible with the results obtained in Figs. [17] and [18] 
In summary, the data seem to indicate a genuine presence 
of small-scale inhomogeneities. Although a specific, model- 
dependent analysis would be needed to better understand 
their origin, the evidence we have suggests that it is not 
unlikely that these substructures are associated to steep 
gradients (or possibly unresolved clumps) in the extinction 
map. 

3.5. Cloud structure functions 

The structure functions of the extinction map of a molecu- 
lar cloud are defined as 



S P (4>) = (\A K (0)-A K (6 + <I>)\ P ) , 



(25) 



where the average is carried over all positions 6 and all 
directions for <fr (note that for p = 2, the structure func- 
tion S2{4>) is the usual two-point correlation function of the 
extinction map Ak)- 



Starting from Kolmogorov (19411, one of the focuses 
of turbulence theory have been the statistical properties of 
the velocity field. In a wide range of length scales, known as 
"inertial range," the effects of both external forces (which 
are taken to act on very large scales) and viscosity (which 
plays a role on very small scales) are negligible. As a re- 
sult, in this range the energy of large scale flows is merely 
transferred to smaller scales until viscous effects become im- 
portant, a process named "energy cascade." This naturally 
leads to random, isotropic motions (because the imprint of 
the large scale flows is likely to be lost during the energy 
cascade), which can be studied with the help of velocity 
structure functions, de fined similarl y to S v (4>) of Eq. (251. 
In his seminal paper, |Kolmogorov| ( 1941[ ) considered the 
first two structure velocity functions S2 and S3, and showed 
that both are simple power laws of the separation <f>, with 




0.0001 



0.0010 



S 3 (4>) 



Fig. 19. The structure functions for the Ophiuchus com- 
plex. The diamonds show the measured S p (0) as a function 
of S3 (</>), with p increasing from 1 to 20; the dashed lines 
connect diamonds with the same p and different c/>; the solid 
lines are the best exponential fits. The first four values of 
p are marked at the left of the corresponding fit. 
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Fig. 20. The scaling index ratio rj(p) /rj(3) as a function of 
the stru cture function order p for the Ophiuchus cloud [cf. 
Eq. ( 26 1] , with the predictions from three theoretical mod- 
els. The bottom panel shows the ratio between the observed 
scaling index and the one predicted by Boldyrev (2002). 
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Fig. 21. Same as Fig. 



20 



but for the Lupus cloud. In this 



case the She & Leveque ( 1994 1 model seems to able to re- 
produce very well the data. Note that no significant differ- 
ence in the scaling index is observed when the low galactic 
latitude regions of the cloud are discarded in the analysis. 



exponents 2/3 and 1. Since Kolmogorov's theory implic- 
itly assumes that turbulence is statistically self-similar at 
different scales, one can actually extend this result to any 
order p, and show that the structure functions of the veloc- 
ity field must be simple power laws of the angular distance 
parameter <j), i.e. 

S p {4>) cx r (p) , (26) 

where T](p) — p/3 in the simplest turbulent model consid- 
ered by Kolmogorov. However, experiments and numeri- 
cal simulations have shown that, although Eq. (26) applies 
to a variety of turbulent flows at high Reynolds number, 
rj{p) substantially deviates from linearity at higher orders 
p, a phenomenon often referred to as "intermittency." |She| 
& Leveque ( 1994 1 have proposed a model for incompress- 
ible turbu l ence b ased on intermittency, further extended by 
Boldyrev (20021 and t ested with nume rical simulations of 
supersonic Turbulence ( |Boldyrev et~aT|2002| ). 

Although no equivalent models are available for the 
(proje c ted) d ensity of molecular clouds, we follow |Padoan 
et al. ( |2002[ ) and assume that the same scaling law rj(p) 
applies to both the velocity and density fields. We analysed 
thus the structure functions S p ((j)) and scaling law i](p) 
for both the Ophiuchus and Lupus complexes. Figure [19 
shows the observed structure functions up to p = 20 on the 
Ophiuchus complex, together with the best exponential fits 
(which are found to be appropriate in our case). The de- 
pendence of the exponent f]{p) of the best fit on the struc- 



ture function order p is shown in Fig. |20[ together with the 
predictions of the three turbulent models discussed above. 
Note that in this anal ysis we cons i dered the ratio ? ?(p)/t?(3), 



1994 
smal 



should show a universal behaviour also at relatively 
Reynolds numbers. While the Boldyrev model seems 
to fit reasonably well the data for the Taurus molecular 
cloud (Padoan et al. 2002 1, the fit seems to be very poor 
for the Ophiuchus and Lupus complexes (Figs. 20 and 21 ). 



Despite the fact that the interpretation of these results 
can be complicated by many factors, the data analysed here 
seems to indicate that the Ophiuchus and Lupus clouds 
have two intrinsically different structure functions. In addi- 



tion, Figs. 20 andpT]s how that one of th e currently favoured 
turbulent models, the Boldyrev ( 2002 ) model, can not de- 
scribe accurately the large-scale structure of these molecu- 
lar clouds. Interestingly, for the Lupus complex a good de- 
scripti on of the scaling index r atio r](p) /rj(3) is instead given 



by the She fc Leveque| ( 1994 1 model, with relatively small 
deviations over the whole range of p values investigated 
here. Note that, by construction, the structure function in- 
dex r]{p) is left unchanged by a simple linear (affine) trans- 
formation Ak {0) <—> aAK(9)+(3 of the extinction map, and 
thus is insensitive to errors on the zero-point of the extinc- 
tion (control field) and on the reddening law. Similarly, the 
structure function index is independent of the distance of 
the cloud: in other words, if two physically identical cloud 
located at different distances d and d! will have structure 
functions that differ only by a scaling factor: 



s'J</>) = s p (k4>) 



where k — d/d' , we have 



(27) 



</>"' (p) cx S' p {<t>) = S p {k<p) cx fc"(p)^(p) cx 



(28) 



or rj'{p) = r](p). Note, however, that if two clouds located 
at different distances and with the same structure func- 
tion index rj(p) (or even with the same structure functions 
S p ) are analysed jointly in the same field, the deduced 
structure functions will in general exhibit a non-trivial be- 
haviour jand will not be simply described as power laws as 
in Eq. ( 26 ) . Since this paper focuses on two wide molecu- 
lar cloud complexes, this point might potentially affect the 
results and partially be responsible for the poor fits dis- 
played in Figs. [20] and [2l| Hence, we decided to repeat the 
analysis in the core of Ophiuchus (see window marked with 
solid line in Fig. [9]) which is most likely comp osed by a 
single clou d at a well determined distance (e.g. |Lombardi 
|et al.|2008[ ) , but the results obtained are completely consis- 



tent to the ones shown in Fig. 20 This, indirectly, confirms 



that for the purposes of the calculation of the cloud struc- 
ture function, there is no significant difference in using the 
whole Ophiuchus complex or a smaller subset centered on 
its core, or equivalently that most likely the poor fit ob- 
served for this cloud with the turbulent models considered 
here is not due to the overlapping, distinct clouds located at 
different distances. Finally, we stress that the conclusions 
reported here seems to be robust against cuts of the star 
catalogue (see Sect. [2| and to changes in the resolution of 
our maps. 
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Cloud 



Total mass Cloud mass 



Ophiuchus 25 400 M Q 8 300 M© 

Lupus 36 500 M© 18 200 M© 

Pipe 14800 M© 14 200 M© 

Pipe (corrected) 9 200 M© 7 600 M© 

Table 2. The masses of the Ophiuchus, Lupus, and Pipe 
dark complexes. For the Pipe, the "corrected" values re- 
ported in the last line refer to the extinction map cor- 
rected for the 0.15 mag plateau. The two columns refers 
to the whole fields considered here, and to the regions with 
Ak > 0.2 mag only. 



4. Mass estimate 

The cloud mass M can be derived from the Ax extinction 
map using the following simple relation 



M = d 2 n(l K [ A K (6)d 2 e , 
Jn 



(29) 



where d is the cloud distance, fi is the mean molecu- 
lar weight corrected for the helium abundance, (3 k — 
1.67 x 10 22 cm~ 2 mag-^ is the r atio [AT(Hi) + N(H 2 )] /A K 
dSavage fc Mathis|[i~979~| see also |Lilley||1955] |Bohlin et aT 
1978^ , and the integral is evaluated over the whole field 
il. Assuming a standard cloud composition (63% hydro- 
gen, 36% helium, and 1% dust), we find \i — 1.37 and to- 
tal masses M = (25 400 ± 2 500) M for Ophiuchus, and 
M = (36 500 ± 3 600) M for Lupus (cf. Table [2). The er- 
ror in both cases is mainly due to the uncertainty on the 
distance of the cloud (fo r both clouds we used our new dis- 
tance measurements, see Lombardi et al. 2008); in addition, 
for Lupus there is an additional uncertainty (not included 
in our error budget) due to the possible projection of other 
clouds on the line of sight. We also considered the total mass 
of the Pipe nebula from our data. As shown by Fig. [6j the 
Pipe nebula is located in an area of fairly large extinction, 
and as a result it is reasonable to subtract from the ex- 
tinction map of the Pipe a constant value, representing the 
"plateau" where the Pipe nebula is located. Clearly, an ex- 
act, physically meaningful definition of the value cannot be 
provided from the available data, because there is no simple 
way to disentangle the effects of different cloud structures 
located at different distances from near-infrared extinction 
measurements alone. Hence, we decided here to derive the 
value of the "plateau" extinction from the same area used as 
a control field in Paper II. This choice is reasonable, because 
it selects the lowest extinction values on the area around 
the Pipe nebula, and in addition produces results that are 
comparable to the ones provided in Paper II. Finally, ad- 
ditionally we considered only regions above 0.2 mag of K- 
band extinction: this choice allows us to concentrate on the 
real structures present in the cloud complexes and to avoid 
the diffuse, low-density material surrounding them. The re- 
sults obtained for the various combinations discussed here 
are reported in Table [2] 

Figure 22 shows the relationship between the integrated 
mass distribution and the extinction in Ak', for the Pipe 
nebula, we used the "corrected" extinction with the plateau 
subtracted. Note that regions with extinction larger than 
Ak > 1 mag account for less than 1% of the total mass 
in the whole field. Hence, we do not expect any significant 
underestimation in the cloud mass due to unresolved dense 
cores. 
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Fig. 22. The cumulative mass enclosed in isoextinction 
contours for the Ophiuchus, Lupus, and Pipe molecular 
clouds. All plots have been obtained from the extinction 
map shown in Fig. [6j and have thus the same resolution 
limit. Note that the difference between the Pipe nebula cu- 
mulative mass function shown here and the one reported 
in Fig. 27 of Paper II is essentially due to the different 
resolution of the two extinction maps. 



More significantly Fig. 22 also shows a clear difference 
between the relative frequencies of extinction characteriz- 
ing the Ophiuchus, Lupus and Pipe clouds. Ophiuchus has 
a considerably higher fraction of dust at high extinction 
(Ak > 1-0 mag) than either Lupus or the Pipe. This re- 
inforces the result of Sect. |3.5| that shows a significant dif- 
ference in the structural properties of the Ophiuchus and 
Lupus clouds at high order p of the struc ture function. 
The scaling index curve for Ophiuchus (Fig. 20 1 is signifi- 
cantly flatter than the corresponding curve for the Lupus 
cloud (Fig. |21~] ) and the predictions of turbulence theory. 
Inspection ofEq. ( 25 ) indicates that this difference is a di- 
rect result of the Ophiuchus cloud having a larger fraction 
of its material at high extinction compared both to that in 
the Lupus cloud or that predicted by standard turbulence 
theory. 

It is interesting to note in this con text that it has been 
known for some time (e.g. Lada 19921 that star formation 



occurs almost exclusively in dense (n(H2) > 10 4 cm -3 ) 
high extinction (Ay > 10 mag) gas. Thus, it is not surpris- 
ing that the Ophiuchus cloud is a very active star formation 
complex, the site of a relatively rich embedded cluster. In 
contrast, Lupus is characterized by mo re modest levels of 
star formation (e.g. Teixeira et al. 2005 1 and the Pipe cloud 
is noted for its nea rly complete absence of star forming ac- 
tivity (e.g. |Alves et al.||2008l . The closer correspondence 
of the Lupus structu re f unction to the predictions of tur- 
bulence theory (Fig. 211 compared to that of Ophiuchus 
(Fig. 20 1 may indicate that standard turbulence inhibits 
the large scale production of dense gas and suppresses ac- 
tive star and cluster formation. 



5. Conclusions 

The main results of this paper can be summarized as fol- 
lows: 



M. Lombardi et al: 2MASS extinction maps: II. Ophiuchus and Lupus 



15 



— We used approximately 42 million stars from the 
2MASS point source catalog to construct a 1 672 square 
degrees Nicer extinction map of the Ophiuchus and 
Lupus dark nebulae. The map has a resolution of 
3 arcmin and an average 2a detection level of 0.5 vi- 
sual magnitudes. 

— We considered in detail the effect of sub-pixel inhomo- 
geneities, and derived an estimator useful to quantify 
them. We also showed that inhomogeneities play a sig- 
nificant role only in the densest cores with Ak > 6- 
8 mag. 

— We derived the structure functions of both dark clouds 
and compared them with several theoretical models. We 
could not find any reasonable fit of the Ophiuchus data 
with models, while the L upus scaling index ratio is well 
described in terms of the She & Leveque ( 1994 1 turbu- 
lent model. 

Acknowledgements. We thank the anonymous referee for many useful 
comments and suggestions. This research has made use of the 2MASS 
archive, provided by NASA/IPAC Infrared Science Archive, which is 
operated by the Jet Propulsion Laboratory, California Institute of 
Technology, under contract with the National Aeronautics and Space 
Administration. CJL acknowledges support from NASA ORIGINS 
Grant NAG 5-13041. 



Miesch, M. S. & Bally, J. 1994, ApJ, 429, 645 
Motte, F., Andre, P., & Neri, R. 1998, A&A, 336, 150 
Murphy, D. C, Cohen, R., & May, J. 1986, A&A, 167, 234 
Nakajima, Y., Tamura, M., Oasa, Y., & Nakajima, T. 2000, AJ, 119, 
873 

Padoan, P., Cambresy, L., & Langer, W. 2002, ApJ, 580, L57 
Padoan, P., Jones, B. J. T., & Nordlund, A. P. 1997a, ApJ, 474, 730 
Padoan, P., Nordlund, A., & Jones, B. J. T. 1997b, MNRAS, 288, 145 
Passot, T. & Vazquez-Semadeni, E. 1998, Phys. Rev. E, 58, 4501 
Savage, B. D. & Mathis, J. S. 1979, ARA&A, 17, 73 
Scalo, J., Vazquez-Semadeni, E., Chappell, D., & Passot, T. 1998, 

ApJ, 504, 835 
Schwartz, R. D. 1977, ApJS, 35, 161 

She, Z.-S. & Leveque, E. 1994, Physical Review Letters, 72, 336 
Stanke, T., Smith, M. D., Gredel, R., & Khanzadyan, T. 2006, A&A, 
447, 609 

Tachihara, K., Abe, R., Onishi, T., Mizuno, A., & Fukui, Y. 2000, 
PASJ, 52, 1147 

Tachihara, K., Toyoda, S., Onishi, T., et al. 2001, PASJ, 53, 1081 
Teixeira, P. S., Lada, C. J., & Alves, J. F. 2005, ApJ, 629, 276 
Vazquez-Semadeni, E. 1994, ApJ, 423, 681 
Vazquez-Semadeni, E. & Garcia, N. 2001, ApJ, 557, 727 
Wichmann, R., Bastian, U., Krautter, J., Jankovics, I., & Rucinski, 

S. M. 1998, MNRAS, 301, L39+ 
Wilking, B. A. & Lada, C. J. 1983, ApJ, 274, 698 



References 

Alves, J., Lombardi, M., & Lada, C. 2008, in Handbook of low mass 
star formation regions, ed. B. Reipurt, ASP Conference Series 

Barsony, M., Kenyon, S. J., Lada, E. A., & Teuben, P. J. 1997, ApJS, 
112, 109 

Benzi, R., Ciliberto, S., Tripiccione, R., et al. 1993, Phys. Rev. E, 48, 
29 

Bertout, C, Robichon, N., & Arenou, F. 1999, A&A, 352, 574 
Bohlin, R. C, Savage, B. D., & Drake, J. F. 1978, ApJ, 224, 132 
Boldyrev, S. 2002, ApJ, 569, 841 

Boldyrev, S., Nordlund, A., & Padoan, P. 2002, ApJ, 573, 678 
Bontemps, S., Andre, P., Kaas, A. A., et al. 2001, A&A, 372, 173 
Comeron, F., Rieke, G. H., Burrows, A., & Rieke, M. J. 1993, ApJ, 
416, 185 

de Geus, E. J. 1992, A&A, 262, 258 
Dubrulle, B. 1994, Physical Review Letters, 73, 959 
Eadie, W., Drijard, D., James, F., Roos, M., & Sadoulet, B. 1971, 
Statistical Methods in Experimental Physics (Amsterdam New- 
York Oxford: North-Holland Publishing Company) 
Greene, T. P. & Young, E. T. 1992, ApJ, 395, 516 
Hughes, J., Hartigan, P., Krautter, J., & Kelemen, J. 1994, AJ, 108, 
1071 

Indebetouw, R., Mathis, J. S., Babler, B. L., et al. 2005, ApJ, 619, 
931 

Johnstone, D., Wilson, C. D., Moriarty-Schieven, G., et al. 2000, ApJ, 
545, 327 

Kleinmann, S. G., Lysaght, M. G., Pughe, W. L., et al. 1994, 

Experimental Astronomy, 3, 65 
Knude, J. & Hog, E. 1998, A&A, 338, 897 
Knude, J. & Nielsen, A. S. 2001, A&A, 373, 714 
Kolmogorov, A. N. 1941, Dokl. Akad. Nauk SSSR, 30, 301 
Krauuter, J. 1991, in European Southern Observatory Scientific 

Report, Vol. 11, Low Mass Star Formation in Southern Molecular 

Clouds., ed. B. Reipurth, J. Brand, J. G. A. Wouterloot, G. H. 

Herbig, B. Pettersson, R. D. Schwartz, L.-A. Nyman, J. Krautter, 

J. A. Graham, B. A. Wilking, & C. Eiroa, 127 
Lada, C. J., Lada, E. A., Clemens, D. P., & Bally, J. 1994, ApJ, 429, 

694 

Lada, E. A. 1992, ApJ, 393, L25 

Lilley, A. E. 1955, ApJ, 121, 559 

Lombardi, M. 2002, A&A, 395, 733 

Lombardi, M. 2005, A&A, 438, 169 

Lombardi, M. 2008, submitted to A&A 

Lombardi, M. & Alves, J. 2001, A&A, 377, 1023 

Lombardi, M., Alves, J., & Lada, C. J. 2006, A&A, 454, 781 

Lombardi, M., Lada, C. J., & Alves, J. 2008, A&A, 480, 785 



